Topological Defect Densities in Type-I Superconducting Phase Transitions 
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We examine the consequences of a cubic term addition to the mean-field potential of Ginzburg- 
Landau theory to describe first order superconductive phase transitions. Constraints on its existence 
are obtained from experiment, which are used to assess its impact on topological defect creation. 
We find no fundamental changes in either the Kibble-Zurek or Hindmarsh-Rajantie predictions. 
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I. INTRODUCTION 

It is generally believed that the Universe, evolving from 
the initial Big Bang, underwent a series of symmetry- 
breaking phase transitions Jl], || accompanied by the cre- 
ation of topological defects, frustrations of the unbro- 
ken phase within the broken one, induced by continuity 
of the order parameter values. These defects appear as 
magnetic monopoles, cosmic strings, domain walls and 
textures. 

Direct experimental tests of these ideas are unfeasi- 
ble, but transitions described by similar equations occur 
in experimentally accessible condensed matter systems, 
and a new trend has unfolded which compares these two 
systems. This "cosmology in the laboratory" relies on the 
fact that the dynamics of phase transitions lie in univer- 
sality classes and the cosmological ones are hence anal- 
ogous to those of condensed matter. For instance, vor- 
tices created in the superfiuid phase transitions of 4 He 
and 3 He have been studied experimentally (see e.g. Ref. 
H |J] for extensive discussions) following an earlier sug- 
gestion in which common features with cosmic strings 
have been noted Q. Similarities between cosmological 
phase transitions and the isotropic-nematic phase transi- 
tion in liquid crystals were studied in Refs. || |7j. Anal- 
ogy with the thermodynamics and transitions in polymer 
chains was drawn in Ref. |^|. 

The case of superconductors is of particular interest 
as the associated phase transition involves a local gauge 
symmetry-breaking process. 

In superconductors, cosmic strings manifest them- 
selves as flux tubes or vortices. Experiments aimed to 
observe defect densities in high-T c materials [g| have lead 
to contradictory results with respect to the density pre- 
dicted by the Kibble-Zurek (K-Z) mechanism This 
prediction is however accurate only for global gauge sym- 
metry breaking, a situation where the geodesic rule for 
phase angle summation is valid. A local gauge treat- 
ment by Hindmarsh-Rajantie (H-R) identifies a new 
mechanism for defect generation, which leads to a pre- 
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diction well below the first Carmi-Polturak experiment 
sensitivity, although the prediction for the second one is 
in reasonable agreement with observation. 

The above experiments were both conducted in type-II 
materials, which exhibit a second order phase transition. 
The question naturally arises as to the extent of changes 
in the defect density predictions for type-I superconduc- 
tors. This is the motivation for our work. 

Distinction between type-I and type-II superconduc- 
tors is traditionally made through the Ginzburg-Landau 
(G-L) parameter k = A/£, the ratio between the mag- 
netic field penetration length, A, and the order parame- 
ter (scalar field) coherence length, £. These characteristic 
length scales are obtained, in the presence of a gauge field 
A, from the free energy density 



F($) = 



2m„ 



iW$ A$ 

c 



(V x A) 



(1) 

where /2 is the magnetic moment of the specimen, m e is 
the electron mass and $ is the order parameter. The G-L 
potential is usually written as 



(2) 



where a is assumed to depend linearly on the tempera- 
ture, a = a (t — 1), t = T/T c , a' and (3 are constants, 
and T c is the critical temperature. Thus one obtains 



A = 



I m e c 2 (3 
Aire 2 \a\ 



(3) 



and 
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FIG. 1: Characteristic Potential curves. 



At T = the coherence length is given by £o = 
%/ \f2m e a' , with k ~ y 7 /?- For k > l/y/2, the transition 
is second order, and £o is typically less than ~ 0.04 \im\ 
for k < l/y/2, the transition is first order, with £o typi- 
cally greater than ~ 0.08 fxm. In general, both are second 
order for H = 0. First order transitions arise from the 
external field term in Eq. ((fl) in the event that a charac- 
teristic sample dimension is greater than A. 

In thermal field theory (TFT) a first order phase tran- 
sition arises from consideration of 1-loop radiative cor- 
rections to a potential of the form of Eq. (^) , which in- 
troduce a barrier between minima through a cubic scalar 
field term, as 



TABLE I: Critical properties of Sn and Al 



Material 


Tc (K) 


F c (0) (G) 


Co (fJ-m) 


A (nm) 


Sn 


3.7 


309 


0.23 


34 


Al 


1.2 


105 


1.6 


16 



II. TEMPERATURE SENSITIVITY BOUND 

A first-order superconductive phase transition mani- 
fests supercritical fields under variation of the tempera- 
ture, as shown in Figure [l]. The superheating curve is 
given by the condition = ^£ = 0, for a certain value 
of $ ^ 0. This yelds a = 9j 2 /16(3. 

In contrast, the supercooling curve is given by the con- 
dition d J^f ^ = for $ = 0, corresponding to a = 0. 
The unobservable critical curve is given by the condition 
V(0) = V(® c ) and = 0, where $ c is the non- 

vanishing minimum of the potential. This corresponds 
to a = 7 2 /2/3. 

Introducing the linear dependencies a = a'(t — 1) and 
7(t) = S t, we obtain for the superheating curve 



9 <5 2 2 

a'Ct - 1) = t 2 

y ' 16 /3 



(7) 



V($) = a$ 2 -7|$|' 
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and hence 



(5) 



where j{T) = ( v / 2/4 7 r)e 3 T O. As in the previous case, 
(3 is constant and a = a!{t— lj is linear with temperature. 

A similar — 7|$| 3 term however, arises in consid- 
erations of gauge field fluctuations in the normal-to- 
superconductor phase transition n3[, with 



7 = Smo^V^Mo^c 
he 



(6) 



resulting in a first order phase transition for all values of 
k. Crossovers between first and second order transitions 
arise in considerations of thermal fluctuations |D| , as also 
in nonlocal BCS treatments 0]. 

With this is mind, we adopt a potential of the form of 
Eq. (^|) and explore constraints on j(T) from experiment. 
The results are compared with TFT 1-loop radiative cor- 
rections. Comparison with the results of Ref. |l2), which 
deals with temperatures close to T c , is attained in the 
limit t — > 1. 

These results are then used to analyze the impact of the 
cubic term on the K-Z and the H-R predictions for type-I 
superconductors, namely on the models themselves and 
on a possible nucleation suppression due to the slowing 
down of the transition induced by the potential barrier. 
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The superheating curve converges to a zero-field shift in 
temperature from T c by (95 2 /16a' (3)T C . This shift is un- 
detected at the current experimental temperature sensi- 
tivity of At exp ~ 10~ 3 p5|. Therefore, a bound for the 
slope of 7 is 



9 6 2 

— —<At e 
16 a'p 



(9) 



The supercooling transition occurs at t = 1 , the critical 
temperature, since it is determined solely by a — (this 
neglects a small correction to a, as given by Ref. |l2|). 



III. SUPERHEATING PERTURBATION 
BOUND 

Measurements of the supercritical fields have com- 
monly been performed on microspheres of type-I mate- 
rials as a means of determining n. Table | indicates the 
critical properties of two, Sn and Al. The existence of 
a cubic term should also be manifest in the presence of 
a magnetic field. To assess its influence, we repeat Ref. 
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p2[ calculations, including the cubic term in the poten- 
tial. For a small superconducting sphere of radius a, the 
magnetic moment is given by p2[ 



a$ A 



3A 2 



a 2 $2 



where $ = ^/^oo and ^ = m e c 2 /47re 2 A 2 . 

After a somewhat lengthy computation (see Ap- 
pendix), the reduced superheating field, h s h = H s h/H c , 
is given by 



1G 



(11) 



where 7g = Z"f /2y/\a\P is a dimensionless parameter, 
and /i^ is the "unperturbed" (7 = 0) superheating field. 

Generally, such measurements have been obtained with 
colloids, and the size distributions of the microspheres 
Jl6[ |i~7| , [l9| pQ, 21 provides a statistical error which ren- 
ders a direct fit of 7g from h s h(t) data unfeasible. Al- 
though the measurement reported in Ref. was con- 
ducted using single microspheres, non-local and impurity 
effects lead to a large theoretical uncertainty in 



h° sh (t) 



^(0)(l-t 2 ) 



(12) 



where /ic(0) = H c (0)/H c (t), which is itself an approxi- 
mation valid only close to T c |23) . 

Since the supercooling field implies the evaluation of a 
second derivative at the origin, it can be easily seen that 
the presence of a cubic term has no effect: 



— ( 7G $ 3 ) = 6 7G $ 







(13) 



The effect of 7 (through j G ) on the superheating field 
must be small, otherwise it would have been already de- 
tected; therefore, we must have 7g < 1, which is not 
valid for relative temperatures in the range 



9 5 2 

l -A—p <t<l • 
For this interval to be vanishingly small, 



(14) 



1a>/3 



< 1 



(15) 



This is a weaker bound than the one of Eq. (|9|). 

For Al and Sn with maximum critical fields of or- 
der 10 2 G, the shift between h 8 h and h° sh is less than 
10 -2 G, well below the sensitivity of measurements 
|L6|, |17|, fjl| [0| |0[ |l|. For this reason, we simply drop 



TABLE II: Derived quantities and bounds for 8 



Material 


Sn 


Al 


a' (J) 


1.15 x 1(T 25 


2.38 x 1(T 27 


P (J.m 3 ) 


4.72 x 1(T 54 


2.16 x 1(T 56 


a' (eV 2 ) 


3.61 x KT 1 


7.45 x 10~ 3 


P 


9.45 x 10" 4 


4.32 x 10" 6 


bound 


t sh shift 


S{eV) 


t sh shift 


S(eV) 




0.25 


1.23 x 1(T 2 


0.25 


1.20 x 10~ 4 




10~ 3 


7.78 x 10~ 4 


10 -3 


7.57 x 10~ 6 


Ref. [12] 


5.19 x 10~ 9 


1.77 x 1(T 6 


2.32 x 10 -6 


3.64 x 10~ 7 


TFT 


2.92 x 10" 9 


1.33 x 10" 6 


3.24 x 10" 6 


4.31 x 10~ 7 



the bound of Eq. (|15| ) an d consider only the one of Eq. 
(Q) . Conversely, a breakdown of the perturbation expan- 
sion of the superheating reduced field would imply a su- 
perheating temperature shifted to t s h = 1.6. Also notice 
that no "spikes" should be seen in the H — T superheat- 
ing curve for values of t in the "exclusion" interval as the 
field values are quite small there. 

Table |0| provides a comparison of the bounds on 5 with 
the prediction of Ref. jl2) . The analogy between cosmol- 
ogy and condensed matter prompts for comparison with 
the TFT cubic term also. To do this, we compute the 
associated slope of j(t) from j(T) — (\/2/47r)e 3 T, ob- 
taining S = (v / 2/47r)e 3 T' c . Obviously, although the pre- 
diction is material independent, its formulation in terms 
of a reduced temperature is not. 

Note that the dimensionality of 7 here is changed with 
respect to the free energy potential of Eq. (||) , through a 
convenient m e factor - this is because the dimension of 
the scalar field in G-L theory is [<I> 2 ] = L~ 3 , its square 
representing a density, while in field theory [$] = L . 
The electron mass determines the conversion as it is ab- 
sent from the kinetic term of the Lagrangean density of 
field theory, d^Qd^Q, while present in the corresponding 
condensed matter free energy term, (?i 2 /2m e )V 2 $ (or, 
equivalently, in the coherence length: £p T = 1/a 1 vs. 



e cm = h 2 /2m e a'). 

Table O also includes the quantities a' and j3, both in 
SI and natural units. As explained above, conversion is 
not direct, but achieved through the multiplicative factor 
m e . 

The results in Table || include the cubic term pre- 
dicted by Ref. [GJJ. In the absence of an applied mag- 
netic field, each momentum-fluctuation of the gauge field 
A has an expectation value given by the equipartition 
theorem. When suitably integrated over the momentum 
space (with a cutoff A of the order of ^ ), 



(A 2 )* = 4^AT c -8n — ^W^T c \<i>\ . (16) 

TT TlC 

Since A 2 couples to $ 2 in Eq. (|l]), this translates into 



4 



an unimportant correction to the scalar field mass, plus 
a negative cubic term, given by — 8/x (e/Tic) ^J-k[IqT c | $ | 3 . 
This term implies a shift in the superheating temperature 
(at zero field), of 



A T = 7.25 x lQ- r2 T*H c (0) 2 $ 



(17) 



with H c (0) in Gauss and £o m M m> 

This shift lies beyond experimental accessibility, since 
it requires a temperature sensitivity of 10~ 6 K (for Al; 
10~ 9 K for Sn). However, such an experiment might be 
performed with Al, using state of the art relative tem- 
perature measurement techniques. 

Surprisingly, for both materials the slopes of 7 pre- 
dicted by TFT and Ref. Q have similar magnitudes, 
~ lO~ 7 eV. This is an indication of the analogous under- 
lying mechanisms behind them: the thermal averaging of 
the gauge field in condensed matter can be thought of as 
equivalent to finite temperature vacuum polarization in 
high energy physics (expressed by the renormalization of 
1-loop Feynmann diagrams). 



IV. TOPOLOGICAL DEFECT FORMATION 

Let us now discuss some possible implications of the 
inclusion of the cubic term in the mean-field potential. 
Since temperature sensitivity measurements constrain 
7g < 10" 2 , we always assume 7 2 <C a(3. 

The K-Z mechanism predicts a density of topolog- 
ical defects (vortices), n ~ (to/ 7 "?)") where t = 
7t?i/16/cbT c is the characteristic time scale, given by the 
Gorkov equation, r q is the quench time, and v is a critical 
exponent. Moreover, it rests upon the assumption that 
there is a single topological defect per £ 2 area and, there- 
fore, one must look for changes induced in this quantity. 
However, since the characteristic scales of the problem 
are obtained via linearization of the G-L equations, close 
to T c and when the order parameter is small, we see no 
changes in this prediction. 

On the other hand, for a thin slab of of width L z , the 
H-R mechanism predicts a defect density of the order n ~ 
(e /2ir)T 1 / 2 LZ 1 ^ 2 i^ 1 , where | ~ 2it/k is the domain size 
immediately after the transition. This quantity is related 
to the highest wavenumber k to fall out of equilibrium 
and is obtained from the adiabaticity relation 



duj(k) 



dt 



\k) 



(18) 



for a given dispersion relation u>(k). In the under- 
damped case, uj(k) = ^Jk 2 + to 2 , with a photon mass 
given by m 2 = 2e 2 |$| 2 = — 2e 2 a//3 . Thus we obtain 

k ~ v/ 'a'e 2 1 ' (3r q , and hence n cx r q . Here, the intro- 
duction of a cubic term in the potential will change the 
photon mass, as the true vacuum shifts to 



$ = 



-37+ \/-16a/3 + 97 2 
4^ 



(19) 



However, since jg = 37/2-v/|a|/3 <C 1, the effect of the 
cubic term is too small to significantly change the H-R 
result. 

Another effect related to metastability concerns the 
non-vanishing probability of the order parameter to 
quantum tunnel from the the symmetric (false) vacuum 
towards the non-symmetric vacuum. Following Refs. 
PH , |25ll , the rate of transition per unit volume and time 
to the true vacua is given, in the thin wall approximation, 
by 



VAt 



2vr T 



3/2 



-Ss/T 



(20) 



where 



Ss(T) 



2ir 1 7 9 (T) 
81 iHT) 



(21) 



is the Euclidean action, and e(T) is the "depth" of the 
true vacuum. 

The thin wall approximation is valid whenever the bar- 
rier's height is much greater than e. This is true when 
7 is comparable to the other parameters, namely when 
7g ~ 1. For eventually smaller values of 7, like those 
predicted by TFT and Ref. [1l2||, the approximation fails. 
In fact, we have shown that the current temperature sen- 
sitivity of 10 -3 -fT only allows values of 7g smaller than 
10~ 2 . Therefore, we must conclude that the barrier's 
height is not comparable to the true vacuum's depth, 
and the field should always tunnel through it (i.e. with 
a probability close to unity). Because of this, there is 
no concern that defects may not have time to nucleate 
within the resolution time of the measuring device, as 
would happen if the potential barrier were high and di- 
minished too slowly. 



V. CONCLUSIONS 

In this work we have examined a possible description of 
a type-I superconductive phase transition by introducing 
a cubic term in the G-L mean-field potential, inspired by 
a gauge field thermal averaging |l2| , and also by analogy 
with TFT. 

Our analysis of the bounds derived from the super- 
heating field and temperature constraints clearly show 
that the contribution of any cubic term is small com- 
pared to other parameters in the G-L potential. Thus 
the following conclusions can be drawn: First, the su- 
perheating temperature shift induced by a cubic term, 
derived either from Ref. 11131 or from TFT, increases 



5 



with decreasing G-L parameter k (At s h(Sn) ~ 10~ 9 ; 
At sh (Al) ~ 10~ 6 ). This suggests that future experiments 
to search for a TFT cubic term should be conducted 
with extreme (k <C 1) type-I materials, for example a- 
tungsten, with T c = 15.4 ±0.5 mK, H c = 1.15 ±0.03 G. 
Similarly, the shift in the supercritical field might be rein- 
vestigated using a DC SQUID, which currently possess a 
sensitivity of 10 _5 ^o/V Hz, or 10 _6 G over a 10 /im grain 
diameter. 

Furthermore, the impact of any cubic term on the de- 
fect density predictions of K-Z @, § or H-R |l§ is neg- 
ligible, with no suppression or slowing down of defect 
production because the potential barrier due to 7 is not 
sufficient to prevent nucleation. 

These considerations suggest that, all else being equal, 
experiments to detect topological defect formation in 
type-I superconductors would observe a reduction of the 
predicted H-R defect densities by 10 — 100 depending on 
choice of material. Recent calculations Q however sug- 
gest that the defect structure formed in type-I materials 
survives significantly longer than in type-II. Given that 
the type-I estimate is of order 10~ 4 seconds, it seems pos- 
sible that the disadvantage in £ might be compensated 
by simple measurability. 



and 



h\ h = (1 + 270*0 - 3*g) j 
Solving for $ we get 



2(1 + 7G $o - *g) = -(1 + 2 7G $ - 3$g) 
which implies that 



/ 3 4 ? 2 



5 + 5 7G 



» 25 lu 5 
Substituting in the first expression, we obtain 



,n 8 . / 3 2 , q . , / 3 2 , , / 3 2 

hl h = -(^- + -7 G ) 3 (1+7 G (V5 + 57 G )-(V5 + 57 G ) 2 ) X 



VI. APPENDIX 



which, to first order in 7c becomes 



Including the presence of a magnetic field with a $ 
dependent magnetic moment, the conditions for the re- 
duced superheating field h s h = H s h/H c are 



(i) minimum 

- 8 - 

h sh - Q 



$g (1 + 7G $o - $0) sirf(io) 



1 + 



sinh(2xo) 
2x 



2sinh 2 (xq) 
^ 



5^15 



'15 



7G 



which is the result Eq. (11) as 



5^15 V A 
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(ii) inflexion point 



$g (1 + 2 7G $ - 3$§) xgsmfe 2 (a;o) 
3 sinh 2 (xo) — Xq coth(xo) — Xq — \xq sinh(2xo) 



where xq = $>oa/\. 

Since the diameter a is much larger than the penetra- 
tion depth A, we can take the limit xq — > 00. The above 
conditions become 
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